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We consider finite-range effects when the scattering length goes to zero near a magnetically con- 
trolled Feshbach resonance. The traditional effective-range expansion is badly behaved at this point 
and we therefore introduce an effective potential that reproduces the full T-matrix. To lowest order 
the effective potential goes as momentum squared times a factor that is well-defined as the scatter- 
ing length goes to zero. The potential turns out to be proportional to the background scattering 
length squared times the background effective range for the resonance. We proceed to estimate the 
applicability and relative importance of this potential for Bose-Einstein condensates and for two- 
component Fermi gases where the attractive nature of the effective potential can lead to collapse 
above a critical particle number or induce instability toward pairing and superfluidity. For broad 
Feshbach resonances the higher-order effect is completely negligible. However, for narrow resonances 
in tightly confined samples signatures might be experimentally accessible. 



PACS numbers: 67.85.-d,03.75.Hh,03.75.Gg,67.85.Lm 
INTRODUCTION 

Recently there has been extended interesting in weakly 
interacting Bose-Einstein condensates for use as an 
atomic interferometer 1; and also to probe magnetic 
dipolar interactions in condensates [2 . This work was 
based on ^^^K atoms where a broad Feshbach resonance 
exists at a magnetic field strength of Bq = 402. 4G [3] 
which allows a large tunability of the atomic interaction 
in experiments [1]. Similar tunability has also been re- 
ported in a condensate of ^Li • The atomic interaction 
can be reduced by tuning the scattering length, a, to 
zero, also known as zero-crossing. In a Gross-Pitaevskii 
mean-field picture we can thus neglect the usual non- 
linear term proportional to a. The question is then what 
other interactions are relevant. As shown in [5], the mag- 
netic dipole will contribute here. 

In the Gross-Pitaevskii picture we might also ask 
whether higher-order terms in the interaction can con- 
tribute around zero-crossing. Recently it was shown that 
effective-range corrections can in fact influence the sta- 
bility of condensates around zero-crossing [6i (7j. The 
Feshbach resonances used thus far in experiments have 
typically been very broad, and as a result the effective 
range, r^, will be small, rendering the higher-order terms 
negligible. However, around narrow resonances this is 
not necessarily the case and finite-range corrections are 
not necessarily negligible. 

The systematic inclusion of finite-range effects through 
derivative terms in zero-range models was begun in the 
study of nuclear matter decades ago [S]. Later on the 
intricacies of the cut-off problems that arise in this re- 
spect was considered by many authors both for the rel- 
ativistic and non-relativistic case (see [H] for discussion 
and references). In the context of cold atoms and Fesh- 
bach resonances, we need to use a two-channel model [10] 
in order to take the lowest order finite-range term into 
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FIG. 1: Scattering length and effective range for the s-wave 
scattering of fermionic ^"K atoms around the Feshbach res- 
onance at Bo = 202.1 G demonstrating the divergence in 
a coupled-channel calculation (symbols) [13] and in a zero- 
range model (full lines). The difference in the zero-range and 
coupled-channel models is caused by the presence of a bound 
state close to threshold in the open channel. 



account. Similar models were already introduced in [TT] 
and denoted resonance models (see f.x. [12] for a compre- 
hensive review of scattering models for ultracold atoms) . 
We note that whereas resonance models treat the closed- 
channel molecular state as a point boson the model of 
[10) treats the molecule more naturally as a composite 
object of two atoms. In the end the parameters of the 
two models turn out to be similarly related to the phys- 
ical parameters of Feshbach resonances. In Fig. ([T]) we 
show calculations of scattering length and effective range 
for the Feshbach resonance at B = 202.1 G in in 
both a coupled-channel model |13j and in the zero-range 
model discussed here. We see the effective range being 
roughly constant at resonance and then start to diverge 
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at zero-crossing. The zero-range model provides a good 
approximation to the full calculations and for many-body 
purposes it is preferable due to its simplicity. Whereas 
the earlier work of ]TT] considered the regime close to the 
resonance, we will be exclusively concerned with zero- 
crossing. To our knowledge the intricacies of this region 
have not been addressed in the literature in the context 
of Feshbach resonances. 

Around zero-crossing the Feshbach model turns out 
to have a badly behaved effective-range expansion. The 
parameters obtained from the effective-range expansion 
should therefore be used with extreme caution as the se- 
ries is divergent at this point. However, as we show in 
this article, the finite-range corrections obtained from the 
full T-matrix at low momenta via an effective potential 
turn out to be the same as one would naively expect 
based on the effective-range expansion. After introduc- 
ing the effective potential we consider its applicability 
and importance in the case of Bose-Einstein condensates 
and for two-component Fermi gases where the attractive 
nature of the effective interaction at zero-crossing could 
lead to collapse above a certain critical particle number 
or to pairing instability and superfluidity. In general, we 
find that tight external confinement is a necessary condi- 
tion for the higher-order effects to dominate the magnetic 
dipole interaction and be experimentally observable. 



TWO-CHANNEL MODEL 

We consider a two-channel s-wave Feshbach model 
with zero-range interactions [10] for which the on-shell 
open-open channel T-matrix as a function of magnetic 
field, is 
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where Afi is the difference between the magnetic mo- 
ments in the open and closed channel, q is the relative 
momentum of the atoms of mass m, abg is the scattering 
length away from the resonance at magnetic field Bq , and 
AB is the width of the resonance. We can compare this 
to the standard vacuum expression for the T-matrix in 
terms of the phase-shift given by 

— gacot d{q) + laq 

where a{B) = abg ^1 — b^^o ) commonly em- 

ployed single-channel models. From Eqs. g and ^ we 
obtain the relation for the phase-shift 
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We now expand the right-hand side in powers of q as is 
usually done in an effective-range expansion. This yields 
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where reo = —2lT?'/{mABAiiabg) is the background value 
of the effective range around the resonance. From Eq. Q 
we can now read off all coefficients in an effective-range 
expansion with their full B-field dependence. For in- 
stance, the effective range is given simply by = 
— l] , which is divergent when a{B) 0. We 
also clearly see that all the other coefRcients are diver- 
gent in that limit. This is signaled also before doing the 
full expansion in q as the first term in Eq. Q diverges 
at zero-crossing. However, in effective potentials derived 
from the T-matrix these problems are not transparent as 
the lowest order coefficient is proportional to a{B) (see 
Eq. (III). Below we will discuss what kind of constraints 
this introduces on the applicability of the effective-range 
expansion near zero-crossing. We note that similar issues 
were briefly discussed in a different context in |14j where 
an equivalent to Eq. ([6]) below was obtained. 

Let us first consider the low-g limit and compare the 
full T-matrix with the effective-range expansion as zero- 
crossing is approached. Taking the low-g limit of Eq. ([s]) 
at zero-crossing where AB/{B — Bq) = I, we find 
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which diverges as q~^. Therefore the coefficients of the 
expansion in Eq. Q must necessarily diverge in order to 
retain any hope of describing the low-g behavior. Fur- 
thermore, since the expansion is an alternating series 
and therefore slowly converged, we also conclude that 
many terms must be retained for a fair approximation 
at very small but non-zero q. The same conclusion can 
be reached by considering the radius of convergence of 
Eq. Q, which we find by locating the pole in Eq. ([s]) at 
ti^q^ /m — A^{B — Bo — AB). This radius indeed goes 
to zero at zero-crossing. We are thus forced to conclude 
that the effective-range expansion breaks down near zero- 
crossing. 



Effective Potential at Zero-crossing 

Since the effective-range expansion is insufficient we 
consider the full T-matrix in the low-g limit at zero- 
crossing. To lowest order we have 
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Using the expression for r^o, this can be written 



(7) 



Knowing the T-matrix at low q we can now proceed to 
find an effective low-q potential through the Lippmann- 
Schwinger equation 



V = T- TGnV, 



(8) 



where Go ~ {E — Hq + i6)^^ is the free space Green's 
function [T3]. This equation can be solved for T{q, q') oc 
q^ +q''^ (the symmetrized version of the full T-matrix) in 
an explicit cut-off approach [HI I15j and then be expanded 
to order q^ for consistence with the input T-matrix. In 
the long- wavelength limit we can take the cut-off to zero 
[15] and for the on-shcU effective potential we then obtain 
the obvious answer 
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in momentum space. The effective potential in real-space 
is now easily found by canonical substitution (q — ;> — iV) 
and appropriate symmetrization 16J. We have 



V{r) = - 



(10) 



Notice that the Lippmann-Schwinger approach is non- 
perturbative as opposed to the perturbative energy shift 
method [Min]- 



Comparison to Effective-Range Expansion and 
Energy-Shift Method 

Away from zero-crossing one can easily relate the 
effective-range expansion to an effective potential 
through the perturbative energy shift method [9, 16..17j. 
To second order the s-wave effective potential is 
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where the first term is the effective interaction usually 
employed in mean- field theories of cold atoms |15j . In 
terms of a and re, we have 52 — — are/2 [16l [17] 
with the field-dependent a — a{B) and — r^iB). 



At zero-crossing the first term in Eq. (Ill vanishes and 



one might expect the second term to vanish as well. How- 
ever, in the naive effective-range expansion of the two- 
channel model discussed above we saw that diverges 
as a^^ and we therefore have 
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In particular, if we for a moment ignore terms in the 
effective-range expansion, we recover exactly the same 



effective potential as in Eq. ( 10 1 at zero-crossing. The fi- 



nite limiting result in Eq. ( 12 ) shows that the potential in 
Eq. (Ill is well-defined as a — 0, provided that appro- 



priate regularization and renormalization is performed. 
Eq. (11) thus applies equally well at resonance (a — >■ 00) 



where the gradient terms are small and at zero-crossing 
where the lowest order delta function term is unimpor- 
tant. It is thus a well-defined effective potential over the 
entire range of a Feshbach resonance. 

We therefore see that even though the effective-range 
expansion has divergent coefficients at zero-crossing, the 
lowest order does in fact give the same effective potential 
as the full T-matrix if we apply it naively. The effective- 
range expansion should thus be viewed as an asymptotic 
series. However, we cannot use the effective-range expan- 
sion to estimate the validity of the second order effective 
potential since the radius of convergence goes to zero at 
zero-crossing as discussed above. 



RELATION TO EXPERIMENTS 

Above we only retained terms of order q^ in the full 
T-matrix. We now estimate the energy regime in which 
this expression is valid. Demanding that the term be 
smaller than the q^ term gives the criterion 
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We relate this condition to recent experiments with 
bosonic condensates of '^^K working around zero-crossing 
[1]. The resonance used there is very broad [AB = 
— 52G) with Gbg = — 29ao and r^o = — 58ao (oq is the 



Bohr radius). The right-hand side of Eq. ( 13 1 is 2.3 • 10 



eV, corresponding to a temperature of about 3 mK. Since 
the experiments are performed at much lower tempera- 
tures the approximation above is certainly valid. How- 
ever, as Ohn and particularly r^o is small, the front factor 
in Eq. ( 10 1 is also small. The relevant scale of compari- 



son is the outer trap parameter b [B] which is typically of 
order l^m, yielding a vanishing ratio lalgVeol/b^ ~ 10"^. 
For broad Feshbach resonances the higher-order interac- 
tions can thus be safely ignored. For very narrow res- 
onances the situation potentially changes as Veo can be 



very large and make the potential in Eq. ( 10 1 important. 
As an example, we consider the narrow resonance in "^^K 
at Bo = 25.85G with AB = 0.47G, Ubg = -33ao, and 
reo = — 5687ao [3] . The right-hand side of Eq. ( 13 1 is now 
2 • 10~^ eV, corresponding to 24 ^K. This is again much 
higher than experimental temperatures. A more care- 
ful argument can be made from the energy per particle 
of the non-condensed cloud. Ignoring the trap, we have 
E/N = 0.nOkBTc{T/Tef/^ {Tc is the critical tempera- 
ture) [in]. For a sample of 3 • 10"* a critical temperature 
of 100 nK was reported in 4J. Using this Tc we find that 
T < QOOnK for Eq. ([Ts]) to hold. Again this is within the 
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experimental regime. The effective potential approach 
should therefore be applicable around zero-crossing for 
narrow resonances. However, even with this narrow res- 
onance we find |a^greo|/fo^ ^ 10"'' and the effect is still 
completely negligible. 

In order to increase the relevance of the higher-order 
term, we now consider some very narrow resonances that 
have been found in ^''Rb. In particular, the resonance at 
Bf) ~ 9.13G [IH] which was recently utilized in nonlinear 
atom interferometry [inj. We have AB = 0.015G, — 
99.8ao, and A^i = 2.00/LtB HHj, which gives reo = -19.8 • 
lO^ao and a ratio \algreo\/b^ = 2.92 • 10-5(l/xm/6)3. A 
trap length of 6 ^ 0.5/im as used in [19] would thus yield 
10~^ and demonstrates that higher-order corrections can 
safely be neglected. For a ratio of 1 we need b ^ 0.03/im 
which is unrealistically small in current traps or optical 
lattices. However, a resonance of width AB — 0.0004G 
is known in the same system at Bq = 406. 2G [21] with 
Qbg = lOOao and = 2.01/ub [20]. In this case we find 
Teo = —7.4 • lO^ao and a much more favorable ratio of 
Wbg^eol/b'^ — 0.001(l/Ltm/&)'^. Here we see that a ratio 
of 1 is achieved already for b ~ 0.1/^m which not far off 
from tight traps or optical lattice dimensions. In terms of 
temperature we still have to be in the ultralow regime of 



T < 30nK according to Eq. ( 13 1 for the latter resonance. 

Consider now a fermionic two-component system 
where s-wave interactions are dominant. Since we have 
< for all Feshbach resonances [2U], the effective 
potential in Eq. ^ is attractive and the system could 
potentially be unstable toward a paired state or become 
unstable to collapse above a critical particle number. 
For simplicity we will use the semi-classical Thomas- 
Fermi approach to describe a gas with equal population 
of the two components and estimate the critical parti- 
cle number. Assuming an isotropic trapping potential 
with length scale b — y/h/muj where w is the trap fre- 
quency, the ground-state density, p{x), can be found by 
minimization and satisfies 
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where p{x) = /cF(a;)/67r^ and a 

imum allowed momentum and chemical potential, /i, is 
found by solving for the turning point of the right-hand 



side of Eq. ( 14 1 which gives 
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We can now compare this k^ax to the value obtained 
from the non-interacting density within the Thomas- 
Fermi approximation at the center of the trap. In terms 
of the number of particles in each component, N, at the 
center of the trap we have fci^(0)& w 1.9067Vi/^ [TS]. By 
equating these two expression we obtain an estimate for 



the critical number of particles, 
evant units, we have 



Nmax- Inserting the rel- 
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where oq is the Bohr radius. We note that the scaling 
Nmax oc can also be obtained by considering the 
point at which the monopole mode becomes unstable. 

Typical numbers for common fermionic species ^Li or 
""^K in the lowest hyperfine states [5D] lead to N^ax ~ 
10^^ for 6=1 /im. This is of course a huge number and 
experiments are well within this limit. Even if one re- 
duced the trap length by a factor of ten and made the pre- 
sumably unrealistic assumption that the particle number 
remains the same we still have N ^ N^ax ■ The reason is 
that the s-wave Feshbach resonances utilized in the two- 
component gases are generally broad in order to study the 
universal regime. If we consider the narrow resonance at 
Bo = 543. 25G in ^Li [22] with AB = O.IG, Ubg = 60ao, 
and Ap = 2.00/iB [20], we have N^ax - 2-10^^{b/lfimf . 
This is somewhat better but we still need b ^ 0.06/im 
to get to an experimentally relevant Nmax ~ 10^. We 
have to conclude that higher-order s-wave interactions 
are highly unlikely to be observable through monopole 
instabilities. In light of this it seems better to consider 
p-wave resonances which are much more narrow in gen- 
eral. However, also here extremely small trap sizes ap- 
pear necessary [^. 

The instability toward Cooper pairing around zero- 
crossing can also be estimated in simple terms. In gen- 
eral the critical temperature is ^ Tp exp(— l/A^olf/j), 
where Nq — mkp{0) /2Tr'^ fi^ is the density of states at the 
Fermi energy in the trap center and U < is a measure 
of the attraction. For the latter we use the effective po- 
tential in momentum space from Eq. ^ and make the 
assumption that q ^ kp(0). Using the expression for 
kp(0) in terms of N above, we find 
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atg J \reo\ 

"'"K this exponent is of 
However, 



For broad resonances in ^Li or 
order 10^ and Tc is thus vanishingly small 
the scaling with trap size can help and if we imagine 
reducing to 6 = 0.1 /im, we find < O.bTp for N = 10^ 
atoms. For the narrow resonance in ^Li discussed above, 
we find that Tc ^ 0.5Tf with N = 10^ can be achieved 
for b ~ 0.5/(m and - O.lTi^ for N = 10^. Thus there 
may be a possibility to reach the pairing instability near 
zero-crossing if high particle numbers can be cooled in 
tight traps and narrow resonances are used. 



Dipole-Dipole Interactions 

The discussion above ignores the dipole-dipole inter- 
action discussed in the introduction which will compete 
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against the higher-order effective potential from the Fesh- 
bach resonance. A simple estimate can be made along the 
lines of the discussion in |15j . The external trapping po- 
tential is the characteristic scale of spatial variations and 
we thus find a ratio, r, of magnetic dipole-dipole, Umd, to 
higher-order s-wave zero-range interaction strength, 1/2, 
which can be written as 



U2 



= 35.7 



1/im 



lOOao 

atg 



lOOOao 



(18) 



For r < 1 the higher-order interaction term will there- 
fore dominate the magnetic dipole term. For the case 
of narrow resonances in ^^Rb discussed above we find 
r - 0.11(5/l^m)2 for the resonance at Bq = 9.13G and 
r - 0.05(6/l/^m)2 for the one at Bq = 406.2G. For 
the narrow resonance in ^Li at Bq = 543. 25G we find 
r ~ 1.4(6/l/im)^. These ratios clearly indicate that mag- 
netic dipole-dipole interactions can be suppressed rela- 
tive to higher-order zero-range terms for narrow Fesh- 
bach resonances and standard trap sizes. This domani- 
nace becomes even stronger for the tight traps needed for 
the realization of the effects discussed above and we thus 
conclude that interference of the magnetic dipole-dipole 
term is not a major concern. 



CONCLUSIONS 

In this article we have discussed the effective potential 
around a Feshbach resonances as the scattering length 
is tuned to zero and finite-range corrections become im- 
portant. We showed that the effective-range expansion 
is badly behaved and the effective potential most be de- 
fined from the T-matrix. We have demonstrated that 
the low momenta effective potential obtained from the 
full T-matrix agrees with one obtained naively from the 
effective-range expansion when the scattering length goes 
to zero. Thus even though the effective-range expan- 
sion has divergent coefficients at zero-crossing the first 
terms of the associated effective potential yield consis- 
tent results. We then estimated the effects of the terms 
on different condensates. Since the effective potential at 
zero-crossing is attractive it may induce various instabili- 
ties which we considered for the case of a two-component 
Fermi gas under harmonic confinement. 

For the broad Feshbach resonances used in current ex- 
periments the effective potential discussed here are neg- 
ligible and the dipole-dipole interaction dominates com- 
pletely at zero-crossing. However, for narrow resonances 
in very tightly confined systems some of the effects might 



be detectable. The competing dipole interaction is small 
for narrow resonances in tight confinement. However, it is 
conceivable that effects of spherically symmetric higher- 
order terms could be separated from dipolar effects which 
change with system geometry [2 . 
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